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Abstract 

We characterize the price of a European option on several assets for a very risk 
averse seller, in a market with small transaction costs as a solution of a nonlinear 
diffusion equation. This problem turns out to be one of asymptotic analysis of nonlinear 
parabolic PDE, and the interesting feature is the role of a nonlinear PDE eigenvalue 
problem. In particular, we generalize previous work of G. Barles and H. Soner who 
studied this problem for a European call option on a single asset where the associated 
eigenvalue problem involves an ODE with an explicit solution. 
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1 Introduction 

In their celebrated paper [1], F. Black and M. Scholes derived a formula for the fair price of 
a European call option on a single stock in an arbitrage free market. They also presented a 
"replication portfolio" that enabled the issuer of the option to hedge his position upon selling 
the option. The Black-Scholes model presented the first rational method for valuing options, 
and consequently, this model has been used in a large number of industrial applications. 

Aside from the financial implications, interesting mathematics also came out of this work 
as the Black-Scholes formula is a solution of a certain linear, parabolic PDE 



now known as the Black-Scholes equation. The purpose of this paper is to discuss an exten- 
sion of the Black-Scholes model and further connections between non-linear PDE and option 
pricing. 

The Black-Scholes model, being the first of its kind, has various shortcomings. One such 
shortcoming is the assumption that there are no costs for making transactions; in fact, in the 
Black-Scholes model, the issuer of an option is trying to hedge his position at each moment 
of time and thus transaction costs would be ruinous. This fact has been formalized and 
proved rigorously [21]. Another shortcoming of the model, is that it does not account for 
risk preferences of option issuers or purchasers; option prices are the same for buyers and 
sellers and each price is completely determined by known market parameters and the option's 
payoff. 

An alternative model, that addresses the aforementioned modeling issues, was presented 
by Davis, Panas, and Zariphopoulou [7]. This model (which we will call the DPZ model) 
uses the principle of certainty equivalent amount to define option prices and poses the option 
valuation problem as a problem of stochastic control theory. Within the DPZ model, G. 
Barles and H. Soner |2] discovered that in markets with small proportional transaction 
costs ~ i/e, the asking price z"^ of a European option by a very risk averse ~ - seller is 
approximately given by 



as e tends to p]. Here ip is a. solution of a PDE resembling a non-linear version of the 
Black-Scholes equation 




r^jj = 0, 
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and A and the function u arising in the error term for together satisfy the ODE eigenvalue 
problem: for each A G M, find X{A) and x ^ u = u{x; A) satisfying 

max |a - y (A + A\i" + {x + Auf) , \u\ - l| = 0. (1.1) 

Estabhshing the convergence of 2;^ to ip, as e tends to 0, is a problem of asymptotic analysis 
of parabolic PDE cLS Z is cl solution of the PDE 



max 




Given the above convergence result, it is very natural to ask if similar phenomena occurs 
for European options on several assets. That is, in the DPZ model for option pricing on 
several assets, does the large risk aversion, small transaction option price exist? And if so, 
can it be characterized as a solution of a nonlinear Black-Scholes type of PDE? The purpose 
of this paper is to establish that this is indeed the case. The main novelty of this work is 
our treatment of the analog of eigenvalue problem described above. G. Barles and H. Soner 
observed that equation (11. ip has a near explicit solution, and this seems to be far from the 
case in the several asset setting; see Theorem 11.11 below. 



1.1 Statement of results 

We consider solutions z" = z''{t,p, y) of the backwards parabolic equation 



max 

l<j<r 



X |-2;t - ^tr (^{p)aa'd{p) [dIz + -^{D^z - y) ® {DpZ " 2/))) , K\ - V^Pij = 0,0 

(1.2) 

where {t,p,y) e (0,T) x (0, 00)" x W, that satisfy 



Here 



z{T,p,y) = g{p). 

( Pi 



;i.3) 



d{p) := diag(pi,p2, 



,Pr, 



P2 



\ 



\ 



Pn J 



and it is assumed that e and T are positive and that a is a nonsingular n x n matrix. Our 
goal is to understand the behavior of solutions when e tends to 0. 

In analogy with the aforementioned work of G. Barles and H. Soner [2], we shall see that 



z'{t, p, y) ^ i){t, p) + eu[ d{p) 



D^{t,p) - y 



;i-4) 



^ For a,b ^ M", a® b is the n x n matrix with i, jth entry Uibj 
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as e tends to 0, where i/j is a. solution of the non-linear Black-Scholes equation 

Vt + e-^(^-*) A (e^(^-*)d(p)L>2^d(p)) + rp ■ Di; - r^p = 0, {t, p) E (0, T) x (0, oo)'^ 

^ = 9, G {T} X (0,oo)« 

(1.5) 

As in the single asset case, we will also see that the non-linearity A and the function u 
together satisfy the following PDE eigenvalue problem: for each A E S{n) H find X{A) and 
X ^ u = u{x] A) satisfying 

max I A - ^trcra* [A + AD'^uA + (x + ADu) ® (x + ADu)) , \u,c, \ - l| = 0. (1.6) 

Our first theorem is 

Theorem 1.1. (Solution of the eigenvalue problem) For each A E S{n), there is a unique 
A = X{A) such that (11. 6p has a viscosity solution u E C(M") satisfying 

lim^^^^ = l. (1.7) 

Moreover, associated to \{A) is a convex solution u satisfying fll.7p .- when det A ^ u E 
C];"{W) for each < a < 1. 

It follows from Theorem 11.11 that the eigenvalue problem associated to the PDE (II. 6p has 
a well defined solution A : S{n) — )■ M. In order to properly interpret solutions of (II. 5p . we 
will need to know that A is a nondecreasing function with respect to the partial ordering on 
iS(n). In fact, we show more. 

Theorem 1.2. (Properties of \) Let A : S{n) — )> M 6e as described in the statement of 
Theorem Then 
(i) A is nondecreasing, 
(a) A is convex, 

(Hi) for each A E S{n) and each permutation matrix U , 

\{UAU') = X{A), 

and (iv) \_ < \ < A+, where 

X_{A) = sup I in^f^ ^traa^ {A + AD'^(f){x)A + {x + AD(p{x)) O (x + AD(f){x))) : 



E C"(M"), max < 1 

l<i<n 



} 



^S{n) denotes the set of real symmetric, n x n matrices. 
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and 

\+{A) = inf i sup ^tma^ {A + AD^tp{x)A+ {x + ADij{x)) {x + AD^p{x))) : 

i) e C"(R"),liminf ^t^^,^ , > 1 

Furthermore, X-{A) = X{A) = X^{A), provided detA ^ 0. 

Employing the eigenvalue function A and making some natural assumptions, we establish 
the following theorem, which is the main result of this paper. 

Theorem 1.3. (Convergence of solutions) Assume that g G C((0, oo)"'). Then for each 
e, T > there is a viscosity solution z" G C((0,T] x (0, oo)" x R") of (11. 2p that satisfies 
fll.3p . Further suppose that 

< g{p) < L (1.8) 



or 



gjp) 

n 



< g{p) < LY,V^ and lim = L, (1.9) 



for a given constant L. Then, as e tends to 0, converges uniformly on compact subsets of 
(0,T) X (0, oo)" X to a viscosity solution of equation (II. 5p . 

In section 121 we study the eigenvalue problem in detail and prove Theorem ll.il In section 
131 we prove Theorem 11.21 which verifies some important properties of A. Finally in section |H 
we establish Theorem II. 3| which characterizes the large risk aversion, small transaction cost 
option price. Before undertaking this work, we present the mathematical model from which 
the equations derive and perform some formal computations that will guide our intuition for 
analyzing z"^ for e small. 



1.2 Model and formal asymptotics 

The market model. Following the work of [71 |2l [9] , we consider a Brownian motion based 
financial market consisting of n stocks and a money market account (a "bond" ) with interest 
rate r > 0. The stocks are modeled as a stochastic process satisfying the SDK 

n 

(iP^(s) = ^aijP^(s)rfiy^(s), s>0, i = l,...,n 
i=i 

where {W{t),t > 0) is a standard n-dimensional Brownian motion and a is a non-singular 
n X n matrix. We assume each participant in the market assumes a trading strategy which 
is simply a way of purchasing and selling shares of stock and the money market account. 
Furthermore, in this model we assume that participants pay transaction costs that are pro- 
portional to the amount of the underlying stock; the proportionality constant we use is 
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On a time interval [t, T], a trading strategy will be modeled by a pair of vector processes 
(L, M) = {{L^, . . . , I/"), (M^, . . . , M"')). Here represents the cumulative purchases of 

the ith stock and M*(s) represents the cumulative sales of the ith stock at time s G [t, T]; we 
assume L*, are non-decreasing processes, adapted to the filtration generated by W, that 
satisfy L^{t) = M^t) = for z = 1, . . . , n. Associated to a given trading strategy (L, M) is 
a process X, the amount of dollars held in the money market, and processes Y^, the number 
of shares of the ith stock held, for i = 1, . . . ,n. These processes are modeled by the SDE 

'dXis) = TX{s)ds + (-(1 + ^e)P\s)dU{s) + (1 - ^e)P\s)dM\s)) ^ < ^ < ^ 
dY\s) = dD{s) - dM\s) i = l,...,n 

We assume that for a given amount of wealth w G M, a seller of a European option with 
maturity T and payoff g{P(T)) > has the utility 

U,{w) = 1 - e~'"/^ 
In particular, the seller has constant risk aversion 

-U'Jjw) _ 1 

If the seller does not sell the option, his expected utility from final wealth is 

v'^f it, X, y, p) = sup EU,{X{T) + Y{T) ■ P{T))E 

L,M 

If he does sell the option, he will have to payout g{P{T)) at time T, so his expected utility 
from final wealth is 



v'{t, X, y, p) = sup EU,{X{T) + F(T) ■ P{T) - g{P{m. 

L,M 

Note that since Ue is monotone increasing, v"^ < v'^'^ . We define the seller's price as the 
amount which offsets this difference (and makes the seller "indifferent" to selling the option 
or not) 

+ A„y,p) = v^'^{t,x,y,p). 

See [5] for more on this approach to option pricing and for more on the theory of indifference 
pricing. 

As in the single asset case [7j, we have the following proposition. Part (i) follows directly 
from Theorem 2 and Theorem 3 of [7]; part [ii) follows from basic calculus. 

Proposition 1.4. (i) v'^^v'^'-^ are the unique continuous viscosity solutions of the PDE 
max I Vy^ - (1 + Ve)PiVx, -Vy- + (1 - \/e)piVx, Vt + ]-d{p)(Ta^d{p) ■ DIv + rp ■ DpV + rxv^ 1=0, 

l<i<n Z j 

(1.10) 

•^Here, and below, we are assuming that X{t) = x,Y{t) = y and P{t) = p. 
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for {t,y,p) G (0,T) x M" x (0, oo)", that satisfy 

v''{T, x,y,p) = 1 — exp{—{x + y ■ p — g{p)) /e) and v"'^ {T, x,y,p) = 1 — exp{—{x + y ■ p) / e) . 
(u)Define z'',z^'f implicitly via 

v" = U^{x + y ■ p — z'') = 1 — exp(— (x + y ■ p — z^)/e) 
v^'-^ = Ue{x + y ■ p — z^'-f) = 1 — exp(— (x + y ■ p — z'''f)/e) 

Then z""^ z""'^ are viscosity solutions of (11. 2p satisfying the terminal conditions 

z%T, p, y) = g{p) and z^'^{T, p, y) = 0. 

Important reductions, (a) To simplify the presentation, we set r = 0. However, this is 
done without any loss of generality as the function 

vit,x,y,p) :=^;(t,e-'-(^-*)x,2/,e-^(^-*V) (l-H) 

satisfies the PDE ( ll.lOp with r = 0, provided of course that w is a solution of f ll.lOp . 
Moreover, v(T,x,y,p) = v(T,x,y,p). 

(b) The main virtue of working with the exponential utility function is that the value 
functions typically depend on the x variable in a simple way. Notice that (upon setting 
r = 0) 

X*'"(s) = ^ + ^ {-(1 + v^)^(^^) ■ dL{s) + (1 - V~e)P{s) ■ dM{s)} , t<s<T 
and so f = w% w^''^ satisfy 

v{t, x,y,p) = l + e-^l\v{t, 0, p)-l). (1.12) 

This is convenient as it reduces the variable dependence of solutions of (ll.lOp . Moreover, 
using (I1.12p . it is straightforward to check that 

2;^, z'^'^ are independent of x. 

The large risk aversion, small transaction cost limit. Directly from the definition of 

Ae and the definition of z^, z'^'^ , we see 

K^ = z'- z'^f. 

Consequently, in order to understand the limiting option price it suffices to study \im.^^Q+ z'" 
and lim£_^o+ '^'''^ ■ Therefore, the problem of characterizing the limiting option price is reduced 
to that of asymptotic analysis of solutions nonlinear parabolic PDE. 

Below, we give a step-by-step formal derivation of how we arrived at the PDE [equation 
(II. 5p ] for the limit i\} and the PDE [equation (II. 6p ] arising in the eigenvalue problem. These 
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heuristic calculations are arguably the most important part of our work since the techniques 
we later use are founded on these results. These computations are based largely on section 
3.2 of P]. 

Step 1. maxi<j<„{|z^. I — i/epi} < 0, so we expect \im^^Q+ z"" to be independent of y. 
This observation leads to the choice of ansatz 

for e small. Here i/), u and are yet to be determined. Using this ansatz, we formally 
compute 

zl ^ ijjt + eDu{x^) ■ xl 
Dyz' ^ e{DyX'fDu{x^) 
Dpz' ^ Dil} + e{DpX'fDu{x') 

D^z' ^ D^iP + e {{Dpx'yD\{x')Dpx' + D^x' ■ Du{x')) 

where (D'^x' ■ Du{x')).. := x^^^, ■ Du{x'), ij = 1, . . . ,n. 
Step 2. We also observe that since 

e\xl. ■ Du{x') \ ^ < y/Ipi, 

x"^ (and its derivatives) should probably scale at worst like 1 / ^/e. With this assumption, we 
calculate 



-zt - itr (^{p)aa^d{p) (^Iz + -^{DpZ - y) ® {DpZ - y) 

-i)t - ^tiaa^ (^d{p)D^i!d{p) + {y/eDpX'd{p)yD\{x'){y/lDpX'd{p)) + 
d(p)£t-l + {^eDpX'd{p)yDu{x')^ ® (^d{p)^^i^ + {V~eDpX'd{p)yDu{x')^^ 



Step 3. Notice that 



V~eDp (^d{p)^^^i^^ d{p) = d{p)D'^d{p) 



( pi^^ 



This basic observation and the above computations lead us to choose the new "variable" 

:= d{p) — 



and the new "parameter" 

A := d{p)D^'ipd{p). 
We further postulate that there is a function A such that 

V^f + X{A) = 0. 

Step 4. With the above choices and postulate, 

r ^ \{A) - ^traa* [A + AD\{x')A + {x' + ADu{x')) ® {x' + ADu{x'))) 
and also for i = 1, . . . ,n 

\u.A^')\<i. 

Since 

max{r, l^j^J - y/epi} = 0, 

l<i<n 

we will require that u and X{A) satisfy 

max I A - -tro-cr* (A + AD'^uA + (x + ADu) (g) (x + , | - 1 1 = 

l<j<n 12 J 

for X G M". In view of estimates we will later derive on 2;^ (see inequality fl4.10p ). we addi- 
tionally require 

1. u{x) 

|x|^oo l^i^^ \Xi\ 

In summary, we have the following nonlinear eigenvalue problem: 

For A E 5(n), find A G M and m : M satisfying 

maxi<i<„ {A - |trcro-* {A + AD'^uA + (x + ADu) (g) {x + ADu)), - l} = 0, x G 
lim|3.|^ooii(a;)/Xir=i 1^*1 = 1 

If we can solve the above eigenvalue problem uniquely for a nondecreasing function A, we 
have the solution of the PDE (11. 5p as a candidate for the limit of 2;*^ as e — 0'''. We remark 
that the procedure described above is philosophically similar to the formal asymptotics of 
periodic homogenization. In analogy with that framework, A plays the role of the effective 
Hamiltonian, and the eigenvalue problem plays the role of the cell problem [201 [12]. Finally, 
we note that the same heuristic argument shows that z"^'^ satisfies the PDE (II. 5p except with 
il)\t=T = terminal condition. We formally conclude that lim^_i.o+ -z*^'-^ = and in particular, 

lim Ae = lim z^. 

e-S>0+ e-5>0+ 
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Remark 1.5. It is possible to use the relationship (11. lip to deduce that if ip satisfies (11. 5p 
with r = 0, then 

satisfies (II. 5p (for nonzero r). Thus, without loss of generality, we will set r = to simplify 
computations. 

Modeling Remarks, (a) While Theorem 11.31 does not cover every possible payoff function 
g for a European option on several assets, it covers many that arise in practice. Some 
prototypical examples are: 

Example 1.6. (Basket call option) p i— j- {Yli=iPi ~ 
Example 1.7. (Basket of put options) p (-)■ Yll=i (-^ ~ Pi)^ 

Example 1.8. ("Mixed" call option) p ^ {pi - K)^ + {pi + 2p2 - K)^ {n = 2) 
Example 1.9. (Forward contract based on the geometric mean) p ^ K — iXYi=iPi)^^^ ■ 
(b) In view of the asymptotic expansion (ll.4p and the growth condition (II. 7p . we have 



A'(t,p,|/) ~ V^(t,p) + ^^epi\iljpXt,p) 



i=l 



as e — 7- 0''". Thus for small e, is a sum of the limiting option price plus a term that 
naturally resembles a transaction cost. 

2 A nonlinear eigenvalue problem 

In this section, we prove the first part of Theorem 11.11 which is the statement that the 
eigenvalue problem is well posed. Our methods are largely based on the approach given 
in our previous article [18], however we consider this work a considerable extension. First, 
we give a definition that will allow for clear statements to follow. This definition involves 
viscosity solutions of nonlinear elliptic PDE and we refer readers to the standard sources for 
background material on this concept O [151 [Ij- We shall also employ the notation of [6]. 

Definition 2.1. u G f/S'C(]R") is a viscosity suhsolution of (II. 6p with eigenvalue A G M if 
for each xq G M", 

max |a - ^tr(7C7* {A + AD'^^{xq)A + {xq + AD^{xq)) ® {xq + ADip{xo))) , \D(^{xq) \ - i| < 0, 

whenever u — (p has a local maximum at Xq and ip G C^(M"'). v G LSC{W^) is a viscosity 
supersolution of (II. 6p with eigenvalue /i G M if for each &W^, 



max 



A _ Iw* {A + AD^,lj{y^)A + (yo + AD^{y,)) ® {y, + AD^P{y^))) , \Diiy^)\ - l| > 0, 
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whenever v — ip has a local minimum at t/o ^"^^ ^ C^(M"). m G C(]R") is a viscosity solution 
of f ll.6p with eigenvalue A G M if its both a viscosity sub- and supersolution of f ll.6p with 
eigenvalue A. 

2.1 Comparison of eigenvalues 

We start our treatment of the eigenvalue problem by establishing a fundamental comparison 
principle that will allow us to compare eigenvalues associated with sub- and supersolutions 
of fOD. 



Proposition 2.2. Suppose u is a subsolution of fll.6p with eigenvalue A and that v is a 
supersolution of fll.6p wi/i eigenvalue fi. If in addition 

hmsup l^""] < 1 < hminf ^4^, (2.1) 

then \ < fi. 

We will first give a heuristic proof, which will naturally motivate the rigorous argument 
to follow. This will be a common theme in our work. 

Formal proof: Assume that M,f G C^(M"). Fix < r < 1 and set 

w'^{x) = Tu{x) — v{x), X G M"". 
By (12. ip . we have that limi^^i^oo w^{x) = — oo, so there is Xr G M" such that 

w'^{xr) = sup w'^{x). 

Basic calculus gives 

= Dw'^{Xr) = TDu{Xr) — Dv{Xr) 
> D'^W^{X^) = TD^u{Xr) - D^v{Xr) 

Note in particular that 

l^x,(a;r)| = T\u^X^r)\ <r <1, i = l,...,n 
and since w is a supersolution of (11.60 with eigenvalue /i, 

/i - ^trao-* {A + AD'^v{xr)A + (x^ + (g) (x^ + AL't;(x^))) > 0. 
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As M is a subsolution of (11.61) with eigenvalue A, 

rA - < [(r - l)A + A{TD'^u{xr) - D'^v{xr))A 

+T{xr + ADu{xr)) ® (x^ + ADu{xr)) - (x^ + ADv{xr)) ® (x^ + ADv{xr))] 

< ^tro-cr* [(r - I) A + t{x^ + ADu{x^)) <S) (x^ + y4Du(x^)) 

-{Xr + TADu{Xr)) ® (x^ + TADu{Xr))] 

= ^tracr* [(r - 1)(A + x^ (g) x^) + r(l - T)ADu{xr) ® AL'm(x^)] 

< i(r - l)tm(7*A + ir(l - T)\a^ADu{Xr)\'^ 

1 7^ 

< -(r-l)tM+-r(l-r)|a*Ap. 
We conclude by letting r — )• 1~. 

□ 

Proof. We now employ a "doubling the variables" argument to make this above heuristic 
proof rigorous. The main difference with standard arguments is that we are working on the 
entire space M". In the proof below, we will quote a few basic facts about semicontinuous 
functions established in [6]. 
1. Fix < r < 1 and set 

w'^(x, y) = ru{x) — v{y), x, ?/ G M". 

For (5 > 0, we also set 

25' 



^s{x,y) = ^Ix-yl"^, x,ye 



The inequality 

u!^{x,y) - ips{x,y) = t{u{x) - u{y)) - ^\x - y\^ + Tu{y) - v{y) 



implies 



< (^Vn\x - y\ - ^\x - y\^^ 



lim {w'^{x,y) - ips{x,y)} = -oo. 



x-y]"^] +ru{y) -v{y) 



Therefore, w'^ — ips achieves a global maximum at a point {xs^ys) € M"' x R". 

2. According to the Theorem of Sums (Theorem 3.2 in [6]), for each p > 0, there are 
X,Y E S (n) such that 

= {D.^(ps{xs,ys),X) e .f'~^{Tu){xs), 
12 



I^x^_m y^ = {-DyMxs,ys),Y) e f'^v{ys), (2.2) 



and 



X 
-Y 



<A + pA\ (2.3) 



Here 



6 V In 



Note that (Q implies that X <Y. 
3. Set 

1 X5 - y5 

P5 = 7 — 

r 

lat p5 G J^'^n(x5). Also note thf 

solutions) , 



and note that ps E J^'^u{xs)- Also note that as maxi<j<„ | < 1 (in the sense of viscosity 



max \tps • ej| < r < 10 



l<i<n 

Since w is a viscosity super-solution of ( II. 6p with eigenvalue /i, we have 

/i - ^tiaa' {A + AY A + {ys + rAps) ® (ys + rAps)) > 
by (12. 2p . As u is a viscosity sub-solution of (11.60 with eigenvalue A, 

A - ^traa* (A + A{X/t)A + {xs + Ap^) ® (x^ + Ap^)) < 0. 

Therefore, 

rA - /X < i(r - l)traa*A + ^traa*[A(X - Y)A\ + ^r|(T*(x5 + Ap^)^ - h,a\y^ + rAp^)^ 

< ^(r-l)trfr(TM+^r|(T*(x5 + Ap5)|2-i|a*(y5 + rAp5)|2. (2.4) 

4. We now claim that G is bounded for all small enough 5 > 0. If not, then there 
is a subsequence of 5 — > such that [w'^ — ip5){xs,ys) tends to — oo. Indeed 



{w'' - ips){xs,ys) = {tu{xs) -v{xs)) +v{xs) -viys) 



F5 - ys] 



26 

< {tu{xs) - v{xs)) + \/n\xs - ysl 



\X5 - ys] 



26 

TlS 

< {ru{xs) - v{xs)) + Y (2-5) 

which tends to —oo as 5 — > provided lim5_j.o+ \xs\ = +oo. This would be the case for some 
sequence of 5 —t- if a;^ is unbounded. 



'Here {ei, . . . , e„} denotes the standard basis in R". 
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However, 



{w^ - (P5){xs, ys) = max <^ tu{x) - v{y) - 



F - y\ 



> Tu{0) - v{0) 

> -oo (2.6) 
and thus xs lies in a bounded subset of M". Since 

lim ; )■ 

(5^0+ 26 

(by Lemma 3.1 in [6]) and ?/5 is bounded, the sequence {{xs,ys))5>o has a cluster point 
(xt-,Xt-) for a sequence of 5 — 0. Note also that ps is a bounded sequence so we can also 
assume that 

as 5 — 0, for some maxi<j<„ \pi\ < 1. 
Passing to this limit in (12. 4p gives 

11 1 
rX-fi < -{t -l)tTa(r^A + -r\(r\xr + Ap)\^ - -\a\xr + rAp)\'^ 

< ^(r - l)tiaa'A + i(r - l)\a'xr\^ + ^r(l - r)\a'Ap\^ 

1 Ti 

< -(r-l)tra(T*A + -r(l-r)|(T*A|2 

We conclude by letting r — )• 1~. □ 
The following corollary is immediate. 

Corollary 2.3. For each A G S{n), there can be at most one A such that (II. 6p has a solution 
u with eigenvalue A satisfying the growth condition (11.71) . 

Now that we know that there can be at most one soluton of the eigenvalue problem we 
are left to answer the question of whether or not a single solution exists. We shall see that 
this is in fact the case. To approximate the values of a potential eigenvalue, we study the 
PDE 



max 

l<j<r 



X - ^tracT* {A + AD'^uA + {x + ADu) ® {x + ADu)) , |u^.J - l| = 0, a; G M" 

(2.7) 

for 5 > and small, and seek solutions that satisfy growth condition (ll.7p 

,. w(x) 
, lim I I = 1- 

The goal is to show that the above PDE has a unique solution us and that there is a 
sequence of 5 — O'^ such that 5us{Q) A(A). Moreover, we hope that us — ^^(O) converges 
to a solution u of (11.61) . First, we address the question of uniqueness of solutions of (12.71) . As 
this can be handled similar to the comparison principle for eigenvalues, we omit the proof. 
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Proposition 2.4. Suppose u is a suhsolution of (12 .Tp and thatv is a supers olution of i \2.7\i . 
If in addition 

u(x) v(x) 
hmsup < 1 < hmmf 

then u < V. 

Corollary 2.5. For each A G S{n), there can be at most one solution of i \2.7\i satisfying 

To establish existence, we need sub- and supersolutions with the appropriate growth as 
|x| — oo. 

Lemma 2.6. Fix < 5 < 1 and A E S{n). 

(i) There is a constant K = K{A) > such that 



uix] 



±M-K]\!^,.,m" (2,8) 

is a viscosity suhsolution of fl2.7p satisfying the growth condition (11.71) . 
{ii) There is a constant K = K{A) > such that 

^ ^ - ^, > 1 

is a viscosity supersolution of (12. 7p satisfying the growth condition (II. 7p . 
Proof, (i) Choose K > such that 

(^\xi\-K^ <^tma'A + ^{\a'x\- ^/^\a'A\)\ xeW. 
As u is convex and as maxi<j<„ \u^. \ = 1, if {p,X) G J^'+n(xo) then 



max I Pi I < 1 and X > 0. 

l<j<n 



Hence, 



1 f V 1 

6u{xo)- -tracer {A + AXA + {xo + Ap)(^{xo + Ap)) < I ^ |xo ■ ei| - A' J - -tracrM 



■l\a-\xo + Ap)\^ 



< (^|xo-ei| -a: ) -hiaa'A 



1 
2 



2 

(|a*xo|-v^|a*A|)' 



< 0. 
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Thus M is a viscosity subsolution. 
(a) Choose 



K := max I -trcra* (A + + (J„ + A)x (g) (/„ + A)x) : max < 1 1 

12 l<j<n J 

and assume that {p,X) G J^'~u{xo). If |a:o ■ ej| < 1 for all i = 1, . . . , n m is smooth in a 
neighborhood of Xq and 

'^(^o) = f + ^ 
Dn(xo) = Xo=p 



D^u{xo) = In = X 



Therefore, 



6u{xo)-^tvaa^ {A + AX A + {xq + Ap) (xq + Ap)) > K-hraa^ [A + A'^ + (/„ + A)a;o (g) (/„ + A)a;o] > 0, 
which imphes 

1 



max <! 5u{xo) - -tiaa^ {A + AX A + (xq + Ap) (g) (xq + Ap)) , \pi\ - 1 ^ > 0. (2.10) 

l<i<n 12 J 

Now suppose |xo ■ ei| > 1 for some i G {1, . . . ,n}. u G C^(M"), so pi = Uxi{xo) = xq ■ ei/\xo ■ 
ei\ and in particular \pi\ = 1. Thus (I2.10p still holds, and consequently, m is a viscosity 
supersolution. □ 

As the existence of a unique viscosity solution now follows directly from applying Perron's 
method of viscosity sub- and supersolutions (see section 4 of [6], for instance), we omit the 
proof. 

Theorem 2.7. Fix < 6 < 1 and A G S{n). There exists a unique viscosity solution u = us 
of the PDE (12.71) satisfying the growth condition (11.71) . 

2.2 Basic estimates 

With the existence of a unique solution of (12. 7p . our goal is establish some estimates on 
that will help us pass to the limit as 5 — > 0. A fundamental property of us that we deduce 
below is that it is convex. Other important estimates of us will be derived directly from this. 
The method of proof is virtually the same as in [18] (Lemma 3.7) and originates from the 
work of [T9l. 



Proposition 2.8. us is convex. 

Proof. 1. We first assume u G C'^{W^) and for ease of notation, we write u for us- Fix 
< r < 1 and set 

nrf \ fx + y\ u{x) + u{y) 

r (x, y) = Tu[ - ' \ x,ye W\ 
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We aim to bound from above and later send r — t- 1^. 

2. As u grows like ^1^=1 1^1 ^ straightforward to check that there is 

{xt-^Ut-) maximizing . At this point, 



- l^Du{Xr) 



and 



Thus, 



tDu 



Xr + yr 



Du{xr) = Du{yr). 



Also observe that v ^ C'^{xt- + v,yr + v) has a maximum at f = which implies 

+ D'^u{xr) + D'^u{yr) 



Since, 



> tD'u 



U 



Xi 



2 

Xr + yr 



< r < 1 



for z = 1, . . . , 77,, we have 



6u{z) - 2^^*^^* + AD'^u{z)A + {z + ADu{z)) (g) (z + ADu{z))) =0, z = Xr, yr- 
Set Zt = {xr + yr)/'^i Pr = Du{zr) , aud uotlce 

~ ^ 2 J 2 

< -traa* (A + AD'^u{zr)A + (2^ + ADu{zr)) ® {zr + 

-^trcra* (A + + {Xr + ADu{xr)) ® {Xr + ADu{xr))) 

tr(T(T* (A + AD^u{yr)A + {yr + ADu{yr)) ® (yr + ADu{yr))) 



4 

(r-1) 



traa^A + -trco-* 

2 2 



+ J|fr*(^, + ADu{zr))\^ - \\a\xr + ADu{xr))\^ - ^\cr\yr + ADn(2/0)p 



2 

(r-1) 



4' 



< -^tro-cj*yl +Ua\zr + Apr)\^ - ^\a\xr + rAp. 



< 



2 
2 

(r-1) 



4' 



\^ - ha\y, + rAp^)\' 



traa'A + . ^\ \a'xr\^ + laVH + ^^(1 - r)\(T'Apr\^ 



4 



tr(T(T*A+ -r(l - r)|(T*A|2 
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for each x, ?/ G M". Sending r — 1 , we conclude that 

x + y\ u{x) + u{y) 



u 



< 0, x,yGM". 



2 J 2 

3. To make this argument rigorous, we fix < r < 1 and now set 



and for ?7 > 0, set 



Notice that 



27] 



x + y 



u 



2 )] 27] 

X + y\ u{x) + u{y) 



x + y 



x + y 



1 

27] 



x + y 



X + y\ u{x) + u{y) 



+TU 



+TU 



From our arguments in part 1 above, it follows that 

lim [uT — if)r^{x,y,z^ = —(X) 

and, in particular, that there is {xr],yr], z^j) G M" x x M" maximizing w'^ — iprj. Now it is 
possible to argue analoglusly to how we did in Proposition 12.21 to conclude that C^^x^y) < 
0(1 — r), as r — 1~. See also the proof of Lemma 3.7 in [18]. □ 



Aleksandrov's Theorem (Theorem 1, page 242 [H]) now implies the following corollary. 
Corollary 2.9. is twice differentiable at (Lebesgue) almost every point in M."^. 
Since Us is convex and ug < u given by (12. 9p . we expect 

5us - -trao-* [A + AD'^usA + {x + ADus) ® (x + ADus)) < Sus - -tiaa^A - -|cr*(x + ADus)\'' 

n ^ 

< K + ^\xi\ tiaa^A 

i=l 



< 
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for all X large enough and S G (0, 1), where K is the constant in (12.91) . In other words, if 
maxi<j<„ \dxiUs{x)\ < 1, then |x| < C for some C independent of 5 G (0, 1). The appropriate 
statement in terms of jets is given below. 

Corollary 2.10. There is a constant C = C{A) > 0, independent ofO<S<l, such that if 
\x\ > C and p G J^'^us{x), then maxi<j<ri \pi\ > 1. 

Proof. Choose C = C{A) so large that 

1 1 

Susiz) < -tiaa^A + 2(1^*^1 - v^l^l)^ kl > C 

for < (5 < 1. Recall that J^'~us{x) = dus{x) by the convexity of us (see Proposition 4.7 in 
[T])|^ Moreover, (p, 0) G .P'~us{x), and so 

1 1 

max{5us{x) - -traa^A - -\o-\x + Ap)\^, \pi\ - 1} > 0. 

l<i<n 2 2 



As 



11 11 
Sus{x) - -tiaa^A - -|cr*(x + Ap)\^ < 6us{x) - -tiaa^A - -(|o-*x| - ^/n\A\f < 0, 



maxi<i<„, > 1. □ 
Corollary 2.11. There is a constant C = C{A) > 0, independent of0<6<l, such that 

us{x) = min <^ us{y) + V \xi -yi\\ , x G M". (2.11) 
\y\<c [ J 

Proof. Choose C = C{A) such that 

" 1 1 

K + '^\xi\~ -tiaa^A- -{\(T*x\- ^/E\A\f <0 for |x| > C, 

2=1 

where K is the constant appearing in the definition of u in equation (12. 9p . Also set v to be 
the right hand side of (12. lip . As maxi<j<„ I^^XiM^I < 1 

Us < V 

and V = Us for |x| < C. 

It is clear that maxi<i<„ \vxi\ < 1, and it is also straightforward to verify that as us is 
convex, v is convex, as well. Now let {p,X) G J'^''^v{xq). If |xo| < C, the v = us is a. 
neighborhood of Xq and so 

max i 6v{xo) - -trcrcr* {A + AX A + [xq + Ap) (g) (xq + Ap)) , - 1 1 < 0. 

l<i<n I 2 I 



^dv{x) := e R" : v[y) > v{x) + p ■ {y - x), for all y € M"} 
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If l^^ol ^ then by the convexity of v 

1 1 
5v{xq) - -tro-o-* {A + AX A + (xq + Ap) ® (xq + Ap)) < 5v{xq) - -tma^A - \a\xo + ADp)\' 



< 5 ^n(O) + ^ |xo ■ Cil j - itraa*A 

-\{\<y'x,\-V^\A\f 

n ^ 

< K + |xo ■ e^l - -trcrcrM 



2 

< 0, 



2 



1=1 

1 



while we always have maxi<j<„ < 1. Therefore, w is a subsolution of (12. 7p . and conse- 
quently 

V < Us- 

□ 

Towards establishing an important lower bound on ug, we first observe that us has its 
global minimum value at x = 0. 

Proposition 2.12. G ^^^(O), and in particular u^ achieves its minimum value at x = 0. 



Proof. By Theorem 12.71 

us{x) = us{-x), X G M"" 

as X I— i- us{—x) satisfies (12. 7p and (11.71) . If us is differentiable at x = 0, then Dus{0) = 0. By 
convexity, 

us{x) > usiO) + DusiO) ■ X = us{0), x G W. 
In general (not assuming differentiability at x = 0), we write u = us and set 

u'{x) = r]' *u{x) = r]'{y)u{x - y)dy, xeR"" (2.12) 

where r;^ G C°° is the standard mollifier (see Appendix C of [8] for more on moUifiers). 
Clearly G C°°(R"'). Recall rj"^ is radially symmetric, is supported in the ball B^, and 
satisfies / 77^ = 1 for all e > 0. As u is continuous, u'^ ^ u locally uniformly as e — )■ 0"*". 
One checks that u*^ is convex and also that u'^(x) = u^(— x). From our remarks above, we 
conclude u'^(x) > u'^(O) for x G M". Sending e — )■ 0"^, gives n(x) > n(0) for all x G M". □ 

We conclude this subsection by establishing an crucial lower bound on u^] this lower 
bound is key to establishing the existence of an eigenvalue. 
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Corollary 2.13. There is a constant C = C{A) > 0, independent ofO<S<l, such that 

us{x)>us{0)+('^\xi\-c\ ,xeW. 



Proof. By above proposition us{x) > us{0) for all a; G and so the claim follows directly 
from Corollary 12.111 □ 

2.3 Existence 

We assume that A is a fixed symmetric, n x n matrix and will now establish the existence 
of a unique eigenvalue X{A). Proposition 12.21 asserts uniqueness, so all that is left to prove 
is the existence of an eigenvalue. To this end, we will use the estimates we have obtained on 
the sequence of solutions us: 



\us{x) - U5{y)\ < X;r=l - y^\ 

us{ix + y)/2) < {us{x) + n5(y))/2. 



for X, y G M" and < 5 < 1. 
Define 

\s ■■= Sus{0) 



vs{x) := us{x) - us{0) 



Notice that 

1 



-tiaa^A <\k< K 



2 

and vs satisfies 

\V6{.X)\ < YJ'i=i 

\vs{x) - vs{y)\ < X;r=i - Vil 
for X G M". We are now in position to prove the existence of an eigenvalue. 

Lemma 2.14. There is a sequence 6k > tending to as k ^ oo, X{A) G M, and u G C(M") 
with \u{x) — u{y) \ < ^11=1 ~ such that 



\\{A) = limfc_^oo A^fe 

1 vsf. u in locally uniformly ask ^ oo 



(2.13) 



Moreover, u is a convex solution of (11.61) with eigenvalue \{A) that satisfies the growth 
condition (11.71) . 
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Proof. It is immediate that X{A) = \imk_^oo ^s^ some Sk 0, as is bounded. The 
convergence assertion of a subsequence vg^. to some u, locally uniformly in M*^, follows from 
the Arzela-Ascoli theorem and a routine diagonalization argument; it is clear \u{x) — u{y)\ < 
Yl^=i \^i~yi\ s-iid that u is convex. It also follows easily from the convergence assertion and 
the stability properties of viscosity solutions (Lemma 6.1 of [6]) that u satisfies the PDE 

max |a(A) - ^tr(T(T* {A + AD'^uA + (x + ADu) ® (x + ADu)) , \u^,^ \ - l| = 0, x G 



Xi 



(in the sense of viscosity solutions). As \u{x)\ < Yll^=i 
for all X e M", 

lim sup — j — - < 1 . 

|x|-s>oo Z^i=l 

By Corollary I2.13[ for all |x| sufficiently large 

n 

vs{x) = us{x) - us{0) > ^ \xi\ - C, 

i=l 

for some C independent of < 5 < 1. Thus, 

(x) 

liminf — ■ — j- > 1, 

and so u satisfies (11. 7p . □ 
To complete the proof of Theorem 11.11 we establish the following regularity assertion. 

Proposition 2.15. Let u he as described in the statement of Lemma ( ]2.14p . If del A ^ 0, 

then u e Cl"{W) for each < a < 1. 

Proof. Note that m is a viscosity solution of the PDE 

max {Lu - b{x, Du), H{Du)} = 0, x G M", (2.14) 

where 

:= — ^ti a a^AD'^ip A is uniformly elliptic, as det A ^ 
b{x,p) := —X{A) + |tr(T(T*y4 + ||cr*(x + Ap)\'^ depends quadratically on p 
H{p) := maxi<j<„ \pi\ — 1 is convex 

Now let O C M" be a bounded domain with smooth boundary and consider the PDE 

' Lu' + (3,{H{Du')) = b{x,Du'), xeO 

u" = u, X G do ' 
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where (/3e)e>o is a family of functions {f3e)e>o satisfying 



Pe = 0, Z <0 

A > 0, z > , , 

2.15 

> ^ ' 

l/3.W = ^, z>2e 

For each e > 0, we think of f3^ as a type of smoothing of ^ i— ?■ {z/e)~^; for small e, we think 
of as a smooth approximation of the set valued mapping 



/3o(t) 




t < 
t = 



Since the values of Pe{H{Du'^)) can be large when H{Du'^) > and e small, solutions will 
seek to satisfy H{Du^) < and, in this sense, become closer to satisfying equation (12 .Mp . 
In this sense, solutions u" of fl2.14p approximate u. 

Theorem 15.10 in [TB] imples that has a unique classical solution u'^ of equation f l2.14p : 
the crucial hypothesis here is that h grows at most quadratically in p. A minor modification 
of proof of Theorem 1.1 in [17] can be used to show that for each O' dd O and a G (0, 1) 
there is a constant C = C{0', a) such that 

\u''\c^,<^{o') < C, 

and moreover that there is a sequence {ek)km tending to such that u'"'-' — )■ u in C^^'^{0), as 
k ^ oo. The important assumption here is that H is convex. From whence it follows that 
u G C[q"(0) and since O was arbitrary, u G C[q"(M"'). □ 

Corollary 2.16. Let u he as described in the statement of Lemma f l2.14p . Then 

Vt := \ x dW : max \ux-{^)\ < 1 

I l<i<n ' 

is open and bounded. Moreover, u G C°°(fi) 

Proof. The first assertion follows immediately from the previous proposition and Corollary 
12. 9[ since x h- )■ Du{x) is continuous mapping of into itself. The second assertion follows 
from standard elliptic regularity, as u satisfies a semilinear elliptic PDE on Q (see Theorem 

6.17 uni). □ 
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3 Properties of the eigenvalue function 



In view of Theorem II ■![ the solution of the eigenvalue problem defines a function that we shall 
denote A : S{n) — > M. In this section, we prove Theorem II. 2[ which details some important 
properties of this function. Our basic tool will be the comparison principle described in 
Proposition 12.21 We use this property to show that A is a monotone, convex function. 
Moreover, the regularity result of the previous section will be used to establish minmax 
formulae for A. In order to establish these properties, we will make use of the following 
characterizations of A, which follow immediately from the existence and uniqueness of the 
eigenvalue function. 

Proposition 3.1. Let A G S{n) and assume that X{A) is the solution of the eigenvalue 
problem associated with equation (II. 6p . Then 



\{A) = sup "I^A G M : there exists a subsolution u of (11.61) with eigenvalue A, 

Vjix] I 

satisfying limsup — j — - < 1- (3.1) 

|x|-s>oo Z^j=i Fil J 

and 

\{A) = inf{ /i G M : there exists a supersolution v of ([L6D with eigenvalue /i, 

V (x) I 
satisfying liminf — j — r > 1. > . (3.2) 

The above formulae, manifestations of the comparison principle, will be used below to 
establish monotone upper and lower bounds on the eigenvalue that will be crucial to deduce 
other properties. Before pursuing these bounds, we note a basic symmetry of A. The following 
proposition states that X{A) is invariant under permutations of A. 

Proposition 3.2. For any A G S{n) and permutation matrix U , we have 

XiUAU') = X{A). 

Proof. Let A G S{n) and u(-; A) be a solution of (II. 6p with eigenvalue X{A). Direct compu- 
tation has that 

v{x; UAU') := u{U'x;A) 

is also a solution of (11.61) that satisfies (11.71) for any permutation matrix U ; the key obser- 
vation here is that u^;. = Du{U^x) ■ and W permutes the standard basis vectors. By 
Proposition O, we have X{A) = \{UAU^). □ 
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3.1 Monotone upper and lower bounds 

In this subsection, we prove that A is a locally bounded, nondecreasing, convex function and 
therefore it is necessarily continuous. We first show that the function A is bounded above 
and below by monotone functions that are constructed from Ai : M — )■ M, the solution of the 
eigenvalue problem found by G. Barles and H. Soner [2]. Then we show A is convex by an 
elementary argument. It turns out that any convex function that is bounded above by a non- 
decreasing function is necessarily nondecreasing itself, and therefore we will be able to con- 
clude that A is monotone. This implies, in particular, that the PDE ipt + X{d(j))D'^iljd{p)) = 
is backwards parabolic which will be useful to us in the following section. 

Proposition 3.3. There are monotone non- decreasing functions A, A : S{n) — ?• M such that 

X{A) < X{A) < X{A), 

for all A e S{n). 

In dimension n = 1, G. Barles and H. Soner [2] showed that the eigenvalue problem 
associated to the ODE 

maxjA - y(y4 + A2u" + (x + - l| = 0, xgR (3.3) 

has a unique solution A : M — R that is monotone increasing and continuous. Moreover, 
associated to X{A) is a solution u = u{- ; A) & C(M) that satisfies 

lim # = 1. 

|3;|-s>oo \x\ 

Furthermore, when A ^ 0, u{- ; A) G C^(R). We will need the following variant of this result 
for our purposes. 

Lemma 3.4. (solution of the ID eigenvalue problem) 

(i) For each A G M and a > 0, there is a unique X = Xi{A, a) G M such that the ODE 



max |a - y (A + A^u" + {x + Au)"^), \u\ - a| = 0, 



X G R (3.4) 



has a solution u = ui{- ; A,a) G C(R) satisfying 



lini ^ = a. (3.5) 

|a;|-5.oo |x| 

When A^O, ui{- ■,A,a)e C\R). 

(a) The function A i— >■ Xi{A,a) is continuous and monotone non- decreasing for each a > 0. 
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Proof. Let A : M — R be the solution of the eigenvalue problem associated to the ODE (13. 3p . 
with solution u = u{-; A) for each A G M as described above. It is easy to check that 

Ui{x; A,a) := u{ax;a'^A), a; G M 

is a solution of (13. 4p with eigenvalue 

X{a^A) 



Xi{A,a) :-- 



that satisfies (13. 5p . The uniqueness of Ai follows from the same ideas used to prove Propo- 
sition IMI □ 



We shall use Ai to design A and A in Proposition 13.31 Our main tool will be formulae 
(13. ip and (13. 2p . First, however, we will perform a change of variables and rewrite (II. 6p . For 
a given A G S{n), we may write 

a^Aa = PAP* 



where P^P = /„ and 



A 



V 



a2 



\ 



Note that 



and for 



v\z) 



\xa^ Acj = flj, 

i=l 

= u{x) with z = a^x, 



tia^aAD'^uA 



ti{a'Aa)D'^v{a'Aa) 
tr[PAP*P)2^PAP*] 
tr[AP*P)\PA] 
tr[K^P*D'^vP] 

n 

i=l 
n 

^ a^D^vPei ■ Pci 



i=l 
n 



'^Pei,Pe, 



(3.6) 



i=l 
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and 



tiaa\x + ADu) ® {x + ADu) 



a\x + ADu)\'^ 
a^x + cr^ADu\^ 
z + {a^Aa)Dv\'^ 
PP^z + PKP^Dv\^ 
P^z + KP^Dv\^ 

n 

^(Z ■ PCi + 

i=l 



Making a further change of variables 

w{y) = v{z) with y = P'^z, 

we have that if 

max I A - -tiaa* \A + AD^uA + (x + ADu) (g) (x + ADi 
l<i<n 1^ 2 '- 



n)l - 1 ^ =0, xeW 



then 



max 

l<i<r 



The above PDE is closely related to the equation 
which has the separation of variables solution 

n n 

A = ^Ai(ai, 1) and = ^ a^, 1), y 



i=l 



i=l 



(here and below, Ai and Ui are a solution pair of (13.31) with = 1). These computations 
motivate the following lemma. 

Lemma 3.5. Let A G S{n) and assume that a^Aa = PAP* where P*P = In and A is given 
by dSSD- 

(z) Set 

n n 

\:='^\i{ai,l/\a\y/n) and w{y) ■.= '^Ui{yi;ai,l/\a\y/n) 



i=l 



1=1 
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where Ai and Ui are the solutions of (13. 4 p as described in Lemma 3. 4' Then u : x ^ 
w{{aPyx) is a suhsolution of (11. 6p with eigenvalue A. 
(a) Set 

n n 

X ■=^\i{ai,\a~^\y/n) and w{y) := Ui{yi] aj, \a~^\y/n) 
1=1 1=1 

where Ai and Ui are the solutions of (13. 4p as described in Lemma \3.4\ Then u : x ^ 
W{{aPyx) is a super solution of (II. 6p with eigenvalue A satisfying 

liminf^4^>l. (3.7) 

Proof. We prove the case where det A 7^ 0, so that ni(-; aj, ■) G C^(M) for i = 1, . . . , ra. The 
general case then follows by straightforward limiting arguments and the stability of viscosity 
solutions under local uniform convergence. 
{€) By assumption, w is a solution of the equation 



As 

n 

traa* [A + AD'^uA + (x + ADu) ® (x + ADu)] = ^ (a^ + a- uj^^j^^ + (y^ + a^w^j^J^) 

and 

max |m„ I = max \aPDw ■ eA 

l<j<n * l<i<n ~ 

< \aP\\Dw\ 

< |cr||-Dw| 



< lo-lv^ max \w„ I 
~ ' ' i<i<n'~y' 

< 1, 



we have that 



max i A - itrcrcr* \A + AD^mA + (x + ADu) ® (x + ADm)] 

l<i<n 12 

Thus M is a subsolution of (ll.6p with eigenvalue A. 
{ii) By assumption, W is a solution of the equation 



M^l-U<0, xeM''. 



max < A (c 

l<i<n I 2 ^ ^ 

- - I i=l 
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Notice that 

n 

traa* [A + AD\A + (x + ADu) ® {x + ADu)] = ^ (a^ + alwy^y^ + {yi + aiWy^f) 



and 



max 

l<j<n 



= max I aPDw ■ Cj 

l<i<n 



> 



1 \Dw\ 



> 



> 



n\iaPy^\ 
1 \Dw\ 

1 |L)in| 



1 maxi<i<„ 











which of course imphes 



Icr ^\\/n ( max lux-l — 1 I > max \wy \ — \a '^\\/n. 

\l<i<n J l<i<n ^ 



It follows that 



max |a - ^trcrcr* [A + AD'^uA + {x + ADu) (g) (x + ADu)] 



, |u,J-U>0, xeW. 



Since, 



where y = {aP)^x, and 



we have 



uix) ^ Yl'^=iUi{yi;ai,\a ^\^/n) 



i<n \yi\ 



ui{t; ai, \a ^\y/n) _^ ^ 

lil-S'OO \t\ 



liminf , > 1. 



Hence, m is a supersolution of fll.6p with eigenvalue A that satisfies (13 .Tp . 
Corollary 3.6. Define 

X{A) := trXi{a'Aa,l/\a\y/E) 



□ 
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and 

A(A) := tr\i{a'Aa,\a-^\^). 

for A G S{n). Then 

A < A < A. 

Proof. For a'^Aa = PAP* where P*P = J„ and A is given by (13. 6p . we have 

n 

trAi((T*A(T, l/lalVn) = ^ Ai(ai, l/lalVn) 

and 

n 

trAi((T*yla, lo-^^li/n) = ^Ai(ai, Io-^-^Ia/k). 

j=i 

Therefore, this corollary follows directly from the above lemma and formulae (13. ip and 
(O. □ 

Proposition 13.31 is now established by noting that 

A K> trAi ((T* Act, 1/1 cr I v^) and A ^-^ tIXl{a^Aa,\a\^/n) 

are nondecreasing. This follows from the proposition below, which is proved in Appendix [XI 

Proposition 3.7. Let f G C{M.) be monotone non- decreasing. Then S{n) 9 A i— )■ trf{A) is 
monotone non- decreasing with respect to the partial ordering on S{n). 

Now we turn to the regularity properties of A and show A is convex and necessarily 
continuous. As we mentioned above, this fact will be used to show that A is monotone 
nondecreasing. 

Proposition 3.8. A : S{n) M. is convex. 

Proof 1. Let Ai,A2 E S{n) and set A3 := {Ai + A2)/2. We will first show 

, / Ai + A2 
A3<^— . 

where Aj := X{Ai), i = 1,2,3. Let Ui = u{-] A^) and assume Ui G C^(W^). We shall only give a 
formal proof, as we now have sufficient experience making the type of argument given below 
rigorous with standard viscosity solutions methods. Finally, we also assume |crcr* = /„. 
A simple inspection of the reasoning below will convince the reader that this can be done 
without any loss of generality. 

Fix r > and note that the function 
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has a maximum on R" x M". For simplicity, we denote this point by (xi,X2) and suppress 
the r dependence. As (xi,X2) is a maximizer 



tDu2.{x2.) = Dui{xi) = DU2{X2) 



(3.8) 



and 



< 0. 



\D^u,{x,) 

The above matrix inequahty imphes that 

^-D^u,{x,){^, + 6) ■ (ei + 6) < \d^u,{x,)^, ■ ^1 + \D^U2{x2)i2 ■ 6- 



for each ^i,.^2 G I^"- Therefore, Lemma [B. II (proved in Appendix [B]) imphes the inequahty 

tr [^3x^2^3(2:3)^] < ^tr [A^D\^{xi)A^] + ^tr [A2D^U2{x2)A2] . (3.9) 

2. As 

\dxiUi{xi)\ = \d^M2{x2)\ = T\d^M3{xs)\ < T < 1, i = l,...,n 

we have 

Xi - tr [Ai + AiD\i{xi)Ai + {xi + AiDui{xi)) (g) {xi + AiDu,{x,))] =0, i = 1,2. 
Therefore, using the inequahty ( 13. 9p gives 
Ai + A2 



< (r - l)trA3 + t\x3 + AsDu3{x3)r 

^\xi + AiDui{xi)\'^ - ^\x2 + A2Du2ix2)\^ 
(r - l)trA3 + r|x3p - ^\x,\' - ^|x2p + t\A3Du3{xs)\' 

^\AiDui{xi)\^ - ^\A2DU2{X2)\^ + 2rX3 ■ A3/^M3(x3) 



Xi ■ AiDUi{Xi) - X2 ■ A2DU2{X2). 



Basic manipulations are 



T\x3?-\\x,\'-\\x2? = {r-l)\x,\' 

and using the first order conditions (13.81) gives 

i/i n I m2 n ^ m2 ^I/i n ^ m2 ^^ \ 1^-0^3(3^3)1^ + \ A2Du3{x3)\^ 
T\A3Du3[x'i)\ - -\AiDui[xi)\ - -\A2Du2[x2)\ = r(l - r) 



Xi — X2 
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and 

2rx3 ■ A3L'm3(x3) - xi ■ AiDuxixx) - X2 ■ A2Du2{x2) = ( — - ) ( ^ ^ ^ ) ■ rDusixs) 

< 



Al- A2\ ^ , , ^ Xi - X2 ^ 



2 I \ 2 

DU3{X3) 

\ ^ J 

Combining the previous four inequalities lead us to 

rXs < (r-l)trA3 + (r- 1)1x31 +r(l-r) < C(l-r) 

for some universal constant C, as maxj Idx^Us^Xs)] < 1. We conclude by letting r — t- 1". 
3. Finally, we remark that virtually the same steps can be used to show 

X{sA + (1 - s)B) < sX{A) + (1 - s)X{B), 

for A,B& ^{fi) and < s < 1. Therefore, the argument above which shows that A is 
midpoint convex also shows A is convex. □ 

Corollary 3.9. A : S{n) — )■ M continuous. 

Corollary 3.10. A : S{n) — M monotone nondecreasing. 

Proof. It suffices to verify the general assertion that if /, (7 : S{n) — )■ M with 

g nondecreasing 
/ convex, and 
J<9, 

then / is nondecreasing (In our case, / = A and g = X). 
Suppose that Q G df{Ao) 7^ 0; that is, 

Q-iA-Ao) + f{Ao)<f{A), AeS{n). (3.10) 

We claim Q > 0. To see this, let ^ G M" and set 

A:=Ao-t^0^ 

for t > 0. As is nondecreasing, substituting this A in (13.101) gives 

-tQ^ ■ e + f{Ao) < f{Ao - ® < giA - ® < giA). 

Clearly this inequality holds for all t > if and only if QC, ■ > i.e. Q > and thus / is 
nondecreasing. □ 
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3.2 M in- max formulae 



In this subsection, we prove two minmax formulae for A. To this end, we make use of 
formulae (13. ip and (13. 2p and the regularity of solutions of equation (11. 6p . We remark that, 
while these alternative characterizations of A are interesting, we do not use them in the proof 
of Theorem II. 3[ Therefore, this subsection can be omitted without any loss of continuity. 

Proposition 3.11. For A E S{n), set 

X-{A) = sup I mf^ i^mcr* {A + ADV(x)A + (x + AD(f){x)) O (x + v4D0(x))) : 



l<i<n J 



and 



X+{A) = inf i sup -traa^ {A + AD'^4j{x)A+ {x + AD4j{x)) (g) {x + AD4j{x))) : 



ij e C^(R"),liminf ^t^^,^ , > U • 

M^oo 2^.^-^ \ Xi\ J 

Then 

X4A) < X{A) < X+{A). 
Moreover, equality holds in both inequalities when det A ^ 0. 

Proof. (A_ < A) Fix A G S{n), let G and suppose that maxj \4>Xi\ < 1- Now set 

fi'f'iA) := ini^ ^trcicr* {A + AD^(j){x)A + {x + AD(f){x)) ® (x + AD(j){x))) . 

If /i'^(A) = — oo, then fJ,'^{A) < X{A); if fJ,'^{A) > — oo, by the assumptions on (p and the 
definition of fi'^{A) 

max I /i'^(y4) - -tiaa^ (A + AD^(f)A + {x + AD(j)) ® (x + AD(j))) , 

l<i<n 12 / ' 

By dSH), we still have /i'^(A) < X{A). Thus, 

A_(A) = sup /(A) < X{A). 

(A+ > A) Again fix A G S{n). Now let ip E satisfy liminf '?/'(x)/ ^"^-^ |xi| > 1 and 



U < 0. 



set 



T^{A):= sup hraa^ {A + AD^^{x)A+{x + ADi:{x))(^{x + ADil:{x))) 
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If T^{A) = +00, then t^{A) > X{A); if t^{A) < +00, by the assumptions on ip and the 
definition of t'^{A) 

max i^r'^ (A) - itiW {A + AD^ijA +{x + ADtP) ® (x + AD^)) , - l| > 0. 

By we still have t'^{A) > X{A). Hence, 

X+{A) = infr'^(A) > X{A). 

■41 

(Equality) Suppose that detA 7^ and let u = u{-,A) be a convex solution of (11. 6p 
associated to X{A) that satisfies u G Cjq'"(M") for each a G (0,1). We first claim that 
X{A) < X^{A). To see this we mollify u 

:= Tj'^ * u 

(see (I2.12P for more on mollification). It is a basic fact that if / G Cjq^(M") then 

\nx) - f{x)\ < |/|c<^(o)e", xeO,:={xeO: dist(x,90) > e} 

for each bounded domain O C M". 

Using this simple observation of mollifiers of Holder continuous functions on M", the fact 
that Du G Cj°^(M") fl L°°(M'^), and some elementary computations, we find 

r]'*{x^ \a\x + ADu{x))\^) = \a\x + ADu')\'^ + Oie") 

as e — )■ 0, for x belonging to bounded subdomains of M". Therefore, as u solves the PDE 
(11. 6p almost everywhere on 

X{A) - ^tiaa' {A + AD^'A + (x + ADu') (g) (x + ADu')) < 0(e") 

for X belonging to bounded subdomains of M". 

Also notice that as m*^ is convex and Du*^ is uniformly bounded, the function 

W 3 x^ {A + AD^u'A +{x + ADu') (g) (x + ADu')) 

has a minimum in a ball -Bk(O) for some R that is independent of e > 0. Consequently, 

A_(A) > inf -tiaa* (A + AD^u'A + (x + ADu') (g) (x + ADn')) 

= inf -traa* (A + AD^u'A + (x + ADm^ ® (x + ADu')) 
> A(A) + 0(e"). 



34 



Next, we claim that A (A) > A+(A). An important observation for us will be that 
maxi<j<„ \uxi \ is uniformly continuous on M". This is due to 

lim max = 1, 

which in turn follows from the limit (11 .7^ and the fact that u is convex. An immediate 
corollary of this uniform convergence is that Du"^ converves to Du uniformly on M", where 
u" = r]^ * u. 
Set 

M^'^ := (l + 5)u^ 

where 5 > is fixed, and notice that 

max \u''/{x)\ < \ <^ max \ul. {x)\ < -. 

l<i<n ^ " l<i<n " 1 + 5 

As 1/(1 + 5) < 1, there is p = p{5) > so small such that 



Also, for eo = eo(^) > small enough 



1 



max u^, X < max \u^^[x)\ + p< —— + p = 7, 

l<i<n l<«<n ' 1 + 

provided < e < eo and maxi<i<„ l^xf(^)l < 1- Moreover, there is ei = ei(5) such that 
{x G : max |ux,(a;)| < 7} C fie := {a: G : dist(x,9fi) > e} 

l<i<n 

for < e < ei- This inclusion follows as the set {x G : maxi<j<„ |na;.(x)| < 7} is an open 
subset of fi. 

Hence for < e < min{eo, ei}, we have 

max \u^'\x)\ < 1 ^ X G fie- 

l<i<n ' 

In particular, if maxi<j<„ |m|'^''(x)| < 1, then 

X{A) - ^traa* {A + AD'^u'A + (x + ADu') ® (x + ADu')) = 0(e"), 

as X{A) - itraa* {A + AD'^uA + (x + ADu) O (x + ADu)) = 0, a.e. on fi. 

With the above computations and the fact that fact that u G C°°(fi) we have 

X+{A) < sup ^traa* [A + AD\''^A + (x + ADu''^) ® (x + ADu''^)) 

sup -tiaa* {A + AD'^u'A + {x + ADu')®{x + ADu')) +0(6) 

maxi \ul'^(x)\<l 

= A(A) + 0(e°) + 0(5). 
We conclude by first sending e and then 5—7-0^. □ 
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Conjecture 3.12. A_(A) = \{A) = \+{A) for all A e S{n). 

Remark 3.13. The conjecture would follow from additional properties of X± such as mono- 
tonicity, convexity, or continuity. 

4 Convergence 

In this section, we verify Theorem 11.31 which characterizes lim^.,. Q-f clS cL solution of the 
nonlinear diffusion equation 

^t + Hdip)D'i/jd{p)) = 0, it,p) e (0,r) x (0,oo)". (4.1) 

Here, A : S{n) — )■ M is of course the solution of the eigenvalue problem discussed in previous 
sections. The method of proof is relatively standard in the theory of viscosity solutions and 
goes as follows. We show the upper limit 

z{t,P,y)-= limsup z^{t\p\y') 

e-S>0+ 
{f ,p' ,y')-^(t,p,y) 

is a viscosity subsolution of (14.11) and the lower limit 

z{t,p,y) := liminf z^{t\p\y') 

e-5-0+ 
{t',p\y')^(t,p,y) 

is a viscosity supersolution of (14. ID . 

As 1 and z_ agree at time t = T and satisfy natural growth estimates for large values of 
p (see Lemma [4.21) . we will be able to conclude 

z < z. 

Combined with the definitions above, we will also have z = z =: ip and that z"^ ^ ip locally 
uniformly as e — (Remark 6.2 in [6]). First, let us make a basic observation. 

Lemma 4.1. 1 and z are independent ofy. 

Proof. 1. As \Zy^\ < ^Jepi for i = 1, . . . , n in the sense of viscosity solutions, 

n 

\z\t,p,y{)-z\t,p,y2)\<Y,^^m\{Vl-y2)■e^V (t,p) G (0,T)x(0, oo)^ yi,y2eM". (4.2) 

2=1 
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Therefore, for ?/i,i/2 ^ I^"" 

zit,p,yi) -z{t,p,y2) = limsup z'{t',p\y[)- limsup z'it'.p.y^) 

< limsup 

(i'>p'>l'i,2/2)^(*.P.yii2/2) 

< lim sup <^ ^fepi | (yi - Z/a) " | 

e^0+ I 

{t',p',y[,y'2)^{t,p,yi,y2) 

= 0. 

Hence, z is independent of y. 
2. As 

:= limsup {-z^){t\p\y') 

(i',p',j/')^(t,P,3/) 

and — also satisfies f l4.2p . we conclude that z_ is independent of y by the same argument 
given in part 1. □ 

We are finally in position to prove Theorem 11.31 The technique we will use, known as 
the perturbed test function method, is due to L.C. Evans [TOj [TTj and was first applied to this 
framework by G. Barles and H. Soner [2]. One difference with the option pricing problem in 
several assets is that we must work with nonsmooth "correctors" i.e. viscosity solutions u of 
equation (11. 6p . We will employ a smoothing argument to overcome this difficulty. 

Proof, (of Theorem ll.3p 1. We now proceed to show that 1 is subsolution of (14. ip and z_ is 
supersolution of (14. ip (with r = 0). First assume that 1 — (p has a local maximum at some 
point (toiPo) ^ (0, T) X (0, oo)" and G C°°; for definiteness, we suppose that 

(2-0)(t,p) < (^-0)(to,Po), (po,to) e 57(^0, Po)@ 

We must show 

- Mto,Po) - X{d{po)D^^{to,Po)d{po)) < 0. (4.3) 

By adding (t,p) H- T]{\t — toP + \p ~ PoP) to and later sending t] 0^, we may assume 
that {to,po) is a strict local maximum point for ^ — in i?T-(to,Po) and also that 

detD^(j){to,Po) ^0. 

We fix (5 > and set 

^A^{t,p) := (1 + 5)2rf(p)Z}V(i,p)t^(p) 

:= A^(to,Po) 
x^''{t,p,y) :=(l + 5)d(p)Ml|b£ 

4)''\t,p,y) := (j){t,p) + eu (a;^'^(t,p, y); Aq) 



^Here ^^(s, z) := {{t, x) G M x M" : |t - + |x - < r^}. 
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for {t,p,y) G (0,T) x (0,00)"- x M". We are assuming that m is a convex, C[q"(]R") solu- 
tion of f ll.6p with eigenvalue X{Ao) that satisfies fll.7p and that the set = {x G M" : 
maxi<j<„ < 1} is bounded. We established the existence of such a function u in the 

proof of Theorem 11.11 

2. We claim there is a sequence of positive numbers tending to zero and local maxi- 
mizers (tfc^Pfc, £ -Br(^O)Po) x 1^" of -2^'' ~ ^'^'"'^ such that 

as A; —7- 00. We will use the idea presented in appendix of [3] to prove this. 

Let yo G M" be given and select a sequence e/c — )■ and {t'^.p'^, y'^j {io.Po, Ho) as — )■ 00 
such that 

-0^-^)(4p'„yl,) -> (^-0)(io,Po) 
(recall z is independent of the y variable). By estimate fl4.10p or (14. lip . 



lim sup 



< -5 < 



YJ!=iV^Pi\yi 

locally uniformly in (t,p) G (0, T) x (0, 00)". Thus, z'"'^ — cj)'^'"^ has a local maximum at some 

{tk,Pk,yk) e 57(to,Po) X M" 

for all k sufficiently large. 

Recall z'' is a subsolution of the eikonal equation 



max 

l<i<n 



{\Zy^ \ - V^Pt} = 0- 



and also that u G C/„'"(M"). Thus, if u^^^'^'^) ^ 



|w.,(x^-^)| < |(l + 5)«.,(x^-^)| 

1 



1 



ekiPk ■ ei){l + 5)u,X^"'^')\ 



(^k{Pk ■ ei) ^ 



< 1. 



In either case, \uxi{x^'' )\ < 1 for i = 1, ■ ■ ■ and hence 



\x 



<^k,S I 



(1+5) 



d{pk 



Di){tk,Pk) - Vk 



for all k> 1. Therefore, the sequence yk is bounded. 
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Without loss of generality, we assume that (ti:,Pk,yk) — > (ti,Pi,yi), as A; — )■ oo. Notice 
that 

{z-(j)){ti,pi) > limsup(2;'^- - (j)""^)itk,Pk,yk) 
>limsup(z^^-0^'='^)(tl,pl,i/O 

k—^oo 

= - 0)(to,Po)- 

As e B-r{to,po), it must be that (ti,pi) = (to,Po)- 

3. Since 

max < 1 

l<i<n 

and u is smooth on the open set Q, 

X{Ao) - ^tiaa'lAo + AoD\{x'''^)Ao + {x'''^ + AoDu{x"'^^)) ® {x"^'^ + AoDu{x"^^^))] = 0, 
for each k >1. Computing as we did in subsection 11.21 we arrive at 



> -(pT'' - itr (^dipk)aa'd{pk) + " Vk) ® (Dpcj)^"' - Vk) 

= -Mtk,Pk) + o{l) 

^ ^.. tiaa'lAo + Ao/^'m(x^'=''^)Ao + + AoDu{x"^'^)) ® + AoDu{x"^'^))] 
2{1 + dj^ 

= -Mto,Po)- ^^^^^, A(Ao) + o(l) 

= -0t(to, Po) - (Y:^^ ((1 + 5)'c/(po)^ '0(^0, Po)d{po)) + o(l) 

as /c — oo. Therefore, we let A; — )■ oo and then 5 — O'^ to achieve (14.30 . 

4. Now assume that 2; — has a local minimum at some point (to,Po) G (0, T) x (0, 00)" 
and G C°°; for definiteness, we suppose that 

(2-</))(t,p) > U-0)(to,Po), (PO,^o) e Brito,Po). 

We must show 

- Mto,Po) - X{d{po)D''(l){to,po)dipo)) > 0. (4.4) 

By subtracting (t,p) H- rj^t — + b ~ PoP) from and later sending 77 — O"*", we may 
assume that (to7Po) is a strict local minimum point for 2 — in i?T-(to,Po) and also that 

detD20(to,Po) ^0. (4.5) 
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We fix 5 G (0, 1), and set 

^A^{t,p) := {I - 5fd{p)D^(P{t,p)d{p) 
Ao:= A^(to,Po) 

^(j)''''P{t,p,y) := (j){t,p) + euP {x'''{t,p,y)) 

for {t,p,y) G (0,T) x (0, oo)" x R". Here u'^ := rj'' * u is the standard mollification of 
u = u{-; Aq), where m is a convex solution of (11.61) with eigenvalue \{Aq) that satisfies (11.71) 
and u G Ci^o"(M") for any < a < 1. 

5. We claim there is a sequence of positive numbers and local minimizers 

{tk,Pk, Vk) G K{to,Po) X M" of z'^ - (P'^^^'P such that 

{tk,Pk) (^o,Po), 

as A; — 7- oo. We will argue as above. 

Let yo G M" be given and select a sequence — )■ and (^'^,^1., y'^) ~^ i'to,Po, Uo) as — )■ oo 
such that 

(^^. - y'k) ^ U - 0)(to,Po) 

(recall z is independent of the y variable). By estimate (I4.13p . 

lim inf ^ i^-r^ >5>0 (4.6) 

I'fi^"^ 2^i=iV(^Pi\yi\ 

locally uniformly in {t,p) G (0,T) x (0, oo)" and all e sufficiently small. Thus, z'^'' — (p'^'^'^'P 
has a minimum at some 

{tk,Pk,yk) e 57(to,Po) X 

for all sufficiently large. Moreover, it must be that yk is a bounded sequence for if not then 
(14. 6 p implies 

while 

and the right hand side above is bounded from above. 

Without loss of generality, we assume that (tk,Pk,yk) — ^ (ti,Pi,yi), as A; — ?• oo. Notice 
that 

< liminfiz'" -4>''^'''P){tk,Pk,yk) 

k—>-co 

<liminf(^^'=-0^-^'O(^l,Pl,2/D 

fe— >oo 

= U- (j)){to,po). 
As (ti,Pi) G i?T-(to,Po)! it niust be that (ti,Pi) = (to,Po)- 
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6. We have at the point itk,Pk,yk) 

10^-^1 = |(l-5)v/^(pfc-e.)<(x^-^)| 

< (1 - 5)v^(pfc ■ Cj) 

< V^kiPk ■ ei) 

fori = 1, . . . ,n. Since ^''^ is a viscosity solution of ([L2]) and ^'^^ G C2((0, T) x (0, cx))"xM"), 
we compute as in subsection 11.21 to get 



< -4>t{tk,Pk) + 

- — tW(Ao + AoL'V(x^'=''')Ao + (x^*-'' + AoDuP{x'>^^^)) ® (x^^''^ + AoL>M^(x^'=''^))). 
2(1 — o)'^ 

(4.7) 

As in the proof of Proposition 13.111 we have that 

A(Ao) - Ui(J(j\Aq + AqD\pAo + (x + A^DuP) ® {x + AqDuP)) < 0{p'') (4.8) 

for X belonging to a bound subset of M" (for some fixed < a < 1). We also have that as 
is convex and Du'' is uniformly bounded that 

A(Ao) - ^traa*(Ao + AqD'^uPAq + (x + Ao^'^i'') ® (x + AqZ^m^)) < A(Ao) - ^trCT(T*Ao 

- (|a*x| - v^|a*A|)2 



< 



for all X large enough. Consequently, (14. 8 p holds for all x G M". 
In particular, (14. 8 p and (14. 7p together imply 

< -Mh^Pk) - (^o) + 0(1) + 0(p°) 

= -Ut^.P,) - 7T^^ (^) + 0(1) + 0{p-) 



(1 


-5)2 




1 


(1 


-5)2 




1 



We obtain (14. 4p by letting k ^ oo and then 5, p — )■ 0"*". 

7. In order to conclude the proof, we need to argue that z < z. Direct computation 
shows that the function 



{t,p) := z{t,p) -v(^ + YIp^ 
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is a subsolution of (14. ip (with r = 0) for each 77 > 0. By Lemma 14.21 below, we have if g 
satisfies (11.81) then 

y9 < 2 < ^ < L, 

and if g satisfies (11.91) then 

n 

ip<z<z<Lj2Pi- (4-9) 

4 = 1 

Here ip is the Black-Scholes price (that satisfies the PDE (14.121) ) and is given by 

JM" V / (ZTTj / 

When g satisfies (II. 9p . the exphcit formula above with inequality (14. 9 p gives 

L = hm =^ < hm =^ < lim =h < L. 

It is now straightforward to check that when g satisfies either (II. 8p or (II. 9p . z"^ — z_ has a 
maximum at some (to^Po) ^ (0, x [0, 00)". 
If to = T, then 

Z < Z + T] + Y^pi 

Letting 77 — > 0"*^ leads to the desired inequality, z < z. Now suppose < T and, for now, 
that 'z, z are smooth. From calculus, 

■.?(to,Po) = ..(to,Po) 

d{po)D^z^{to,Po)d{po) < d{po)D^z{to,Po)d{po) 
However, these inequalities would imply a contradiction as 

1 

^0 



- -Zt{to,Po) + Zf(to,Po) 

On 



< \{d{po)D^z{to,Po)d{po)) - Xid{po)D\{to,Po)dipo)) 
= X{d{po)D^z^{to,Po)d{po)) - Kd{po)Dh{to,po)d{po)) 

< 0. 

The last inequality above is due to the monotonicity of A. It is now routine to use the ideas 
in Section 8 of [6] to make the same conclusion without assuming smoothness of z, z. □ 



Lemma 4.2. Let z'" he the solution of (ll.2p described in Proposition \1.4 
(i) There is a universal constant C such that 

n 

ip{t,p) <z\t,p,y) <C + Y,V~m\y^'Cl {t,p,y) e (0,r] x (0,oo)" xM", (4.10) 

1=1 
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provided g satisfies (11. 8p or 

n n 

ip{t,p) < z'{t,p,y) < C J]p, + ^v^p,|y,-C|, {t,p,y) e (0,T] x (0,00)" x W, (4.11) 

1=1 i=l 

provided g satisfies (II. 9p . Here (f is the "Black- Scholes" price 

' ^t + \traa'{d{p)D^^d{p)) = 0, {t,p) e (0,T) x (0, 00)" 



^ = 9, {t,p) e {T} X (0,00)- 
(a) For each < r] < T , there is a K = K{r]) such that 



(4.12) 



'{t,p,y) > J2 V~m\y^\ - KTe, it,p,y) e (0,T - r/] x (0, 00)" x M'^ (4.13) 



for alio < e < 1 /4. 

Proof, (i) We prove (14. lip as the proof for (I4.10p is similar. Let C := L be the constant in 
f O]) and set 2/0 := (C, C, . . . , C) G M". In view of inequahty ffLQ]) . 

^''(^,P,l/o)|t=T = 1 -exp{-{x + y-p- g{p))/e) 

= 1 - exp ^- + -^(j>)j /e 

> 1 — exp(— a;/e) 

As (t,p) — )■ v^{t,p, yo) is a supersolution of the backwards parabohc PDE 

-^Pt - itraa* (rf(p)D2^rf(p)) = 

(recall equation (ll.lOp ). it must be that v''(t,p,yo) = 1 — exp {—{x + yo ■ p — z^)/e) > 
1 -exp(-x/e) for (t,p) G (0,r] x (0,cx))". Hence, z%t,p,yo) < yo ■ p = CYJ'i=iPi and 
consequently, 

n 



' (t, p, y) < (t, p, yo) + ^ v^P. 

n n 

< C^Pi + ^ y/epi\yi - C 



i=l 1=1 



as claimed. 

Similarly, we observe that (t, x,p,y) H- U^{x + y ■ p — (f) is a supersolution of (ll.lOp that 
satisfies v''\t=T = U^{x + y ■ p — g)- It follows that w ^ < ?7e(x + 2/ ■ p — V') and thus z'' > ip. 
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(ii) We follow the proof of Lemma 2.2 in [2] closely. Let r] e (0,T] and g G C°°(]R) be a 
nondecreasing function such that g{s) = for s < and g{s) = 1 for s > 77. Also set 

n 

'^it,P,y) := giT -t)^y/epi\yi\ -Ke{T -t) 

1=1 

for a positive constant K to be determined below. To establish the inequality (14. lip , it 
suffices to choose K so that ip is a subsolution of flL2p . 

As < g < 1, we have \ipy. \ < y/epi for i = 1, . . . ,n. Moreover, 

1 f.,. ... V^o. 1 



Q: = -iJt- ^ti \^d{p)aa'd{p) (^D^^ + -{Dpi; - ® {Dp^ - y) 

"1 2 
= -(^K + g^y/epi\yi\ - — yd{p){Dpilj - y)\ 



i=l 
n 



< -eK + g'^Vepi\yi\ - J-ip ^P^dV^lVil - Vif 

i=l ' ' 1=1 

n 1 " 

< -eK + g'Y, - 7rr-T\^Y.P^y^^^ " 



2e c7- , 
1=1 ' ' 1=1 



< 

for K = K{rf) chosen large enough. K can be chosen independent of e as ey^/ {1 — y/e) < 1/4 
for e G (0,1/4). □ 

Acknowledgements: Part of this work was accomplished during my graduate studies and 
I thank my thesis advisor L.C. Evans for his encouragement. I am also indebted to Y. Yu 
for many interesting discussions related to this work. 

A Monotonicity oi A tr/(A) 

In this appendix, we prove Proposition 13.71 which we restate for the reader's convenience. 

Proposition A.l. Let f G C(M) be monotone non- decreasing. Then S{n) 3 A trf{A) is 
monotone non- decreasing with respect to the partial ordering on S{n). 

We first prove a version of this proposition for smooth functions. 

Proposition A. 2. Let f G C°°(]R) and N G S{n). For each A G S{n), we have 

Dtr[f{A)]N := M^ilMzJlZi^ = tr[f{A)N]. (A.l) 
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Proof. A = OAO*, where 0*0 = /„ and A is diagonal; we shall also write A = {0'ij)i<i,j< 
and A = diag(Ai, . . . , A„). With this notation, 

n 

tr/(A) = ^/(A,). 

i=l 

Now Afc = AOck ■ Ocfc and A = ^"j^j^ aijCiC^y Therefore 



- — = Cie^Oek ■ Ock + A- — Oe^ ■ Oe^ + AOeu ■ t: — Oe^ 
oaij ■' oaij oaij 

= e]Oek ■ epck + 2A0ek ■ -^Oe^ 



OikOjk + 2AfcOefc ■ 7^ — Ocfc 



d_ 



d 

= OikOjk + \k^ — lOcfcl^ 

O ikO jk^, 

since 0*0 = J„. 

Let G S{n). We have from the definition of Dii[f {Ay\N in (lA.ip . the chain rule, and 
the above computations that 

DiT[f{A)]N = 

i,j=l 

i,j=i k=i ■> 

n n 

= I]/'(Afc)$^0,fcO,fcA^,,- 

k=l i,j=l 
n 

= Y.f{\k){0'NO)kk 

k=l 

= tr[/'(A)Ar]. 

□ 

Lemma A. 3. Let f G (^""(IR) be non- decreasing. Then S{n) 3 A ^-^ trf{A) is non- 
decreasing. 
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Proof. Let A,B^ ^i^^) with B > A. From the above proposition, 

tr/(S) - tr/(A) = ±tTf{A + t{B-A))dt 

= [ DtTf{A + t{B - A)){B - A)dt 
Jo 

= I ii[f'{A + t{B-A)){B-A)]dt 

Jo 
> 0, 

as the matrix f'{A + t{B - A)) > for all t G [0, 1] and 5 - A > 0. □ 

Proof. (Proposition (13.71) ) Let denote the standard mollifier of / and suppose A,B^ ^{f^) 
with B > A. By the above lemma, trf^^B) > tr/^(yl) for all e > as is non- decreasing. 
Letting e 0+ implies ti f{B) > trf{A). □ 

B An elementary matrix inequality 

Lemma B.l. Let a,b,c E S{n) and b,c nonnegative definite. Further assume that 

a(ei + 6) ■ (6 + 6) < &6 ■ 6 + c6 ■ 6- (B.l) 

for all ^1, ^2 e M". Then for any Ai, G S{n), 

tr[{Ai + A2Ya\ < tr[bAl] + tr[cAl]. 

Proof. 1. First assume that b and c are positive definite. Then for a given ,^1 G M", we may 
select ,^2 such that 



Also note that for these vectors 



(6 + c)ei = c(ei+6) 
(6 + 0)^1 = 6(^1 + 6) 



Therefore, 



66 ■ 6 = + c)-'c{^i + 6) • {b + c)-ic(6 + 6) 
= c(6 + c)-i6(6 + c)-ic(ei + 6) ■ (ei + 6) 
= c^6(6 + c)-2(ei+6)-(6 + 6), 

where we are using the fact that all the matrices being manipulated are symmetric. 
Similary 

c6-6 = &'c(&+c)-2(ei + 6)-(ei + e2), 
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and so 

b^i ■ ei + c6 ■ 6 = {b'cib + c)"' + c'bib + c)-2) (ei + 6) ■ (6 + 6) 
= (c^& + &'c)(6 + c)-2(6+6)-(6+6) 

= 6c(6 + c)-i(ei + 6)-(ei + e2). 

As the vector sum + ^2, where b^i = c^2, runs through all vectors in R"-, we have by 
assumption (IB.ip 

a < 6c(6 + c)"\ 

Also observe that symmetric matrices obey the Cauchy-Schwarz inequality: 

«i"2 < ^ , di,d2 e S[n). 

With these observations, we have 

tr [{Ai + A2fa] = ti[Aja + 2^1^20 + Aja] 

= tilAja + 2Vb-'^cAiVc''^bA2a + Aja] 

< tr[((4 + b-'c)Al + (/„ + c-'b)Al)a] 

< tr[((J„ + b~^c)Al + (J„ + c~^b)Al)bc{b + c)-i] 

= tr[6c(6 + c)-i(4 + b~'c)Al] + tr[6c(6 + c)^i(/„ + c"'b)Al] 
= tT[bAl] + tT[cAl]. 

2. If fe, c are only nonnegative definite, we replace them by 6 + 6 In and c + Sin for 5 > 
and small. We can then repeat the computation above and pass to the limit as 5 — t- 0^ to 
get the desired inequality. □ 
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